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ABSTRACT 

We present a manifestly rotational invariant formula- 
tion of the matrix product method valid for spin chains 
and ladders. We apply it to 2 legged spin ladders with 
spins 1/2, 1 and 3/2 and different magnetic structures 
labelled by the exchange coupling constants, which can 
be ferromagnetic or antiferromagnetic along the legs and 
the rungs of the ladder We compute ground state en- 
ergy densities, correlation lengths and string order pa- 
rameters. We present numerical evidence of the duality 
properties of the 3 different non ferromagnetic spin 1/2 
ladders. We show that the long range topological order 
characteristic of isolated spin 1 chains is broken by the in- 
terchain coupling. The string order correlation function 
decays exponentially with a finite correlation length that 
we compute. A physical picture of the spin 1 ladder is 
given in terms of a collection of resonating spin 1 chains. 
Finally for ladders with spin equal or greater than 3/2 
we define a class of AKLT states whose matrix product 
coefficients are given by 9-j symbols. 



I) INTRODUCTION 

The matrix product method (MPM) is a variational 
approach appropiate to study the ground state and exci- 
tations of a variety of one-dimensional lattice systems in 
Condensed Matter and Statistical Mechanics. The the- 
oretical and experimental interest of these systems has 
grown spectacularly in the last years, due to the discov- 
ery of interesting and unexpected physical properties in 
spin chains and ladders. 

The basic idea behind the MPM is the construction 
of the ground state and excitations of ID or quasi ID 
systems in a recursive way, by relating the states of the 
system with length N to that of length N —1. This sim- 
ple idea has appeared in the past in different places. First 
of all, in the Wilson's real space renormalization group 
the ID-lattice is built up by the addition of a single site 



at every RG step Q|. This procedure is also used in the 
density matrix renormalization group method (DMRG) 
of White |2j . Other source of the MPM is the well known 
AKLT state of the spin 1 chain j|. This is a simple but 
non trivial example of a matrix product state, which has 
motivated various generalizations as the ones of Klumpcr 
et al Q|, Ostlund and Rommer ||, etc. We shall follow 
in this paper the formulation of the MPM due to the lat- 
ter authors, which is based on the analysis of the fixed 
point structure of the DMRG ground state in the ther- 
modynamic limit Q. A closely related approach is that 
of Fannes et al |6j . The MPM offers an alternative for- 
mulation of the DMRG method in the regime where the 
latter reaches a fixed point after many RG iterations [Q. 

Whereas the DMRG is a purely numerical method, the 
MPM offers the possibility of an analytical approach to 
elucidate the actual structure of the ground state (GS) 
and excitations. The MPM is a standard variational 
method which determines the variational parameters by 
minimizing the GS energy. Minimization problems are in 
general harder than diagonalization ones. In this respect 
the MPM is so far less performant than the DMRG. How- 
ever we believe that the analytical insights gained with 
the MPM could be used to boost the numerical precision 
and applications of both the MPM and the DMRG. 

In this paper we apply the MPM to the 2-leg spin 
ladder. Spin ladders with diagonal couplings have been 
studied with the MPM of Klumper et al in ref. [§|. Spin 
ladders were first studied as theoretical labs to test ideas 
concerning the crossover from ID to 2D, with the surpris- 
ing result that this crossover is far from being smooth: 
the even and odd ladders display quite different prop- 
erties converging only when the number of legs goes to 
infinity (for an introduction to the subject see Q). Even 
spin ladders are spin liquids with a finite spin gap and 
finite spin correlation length, while odd spin ladders be- 
long to the same universality class than the spin 1/2 
antiferromagnetic Heisenberg chain, which has no gap 
and the correlations decay algebraically. Another rea- 
son to study ladder systems is that materials actually 
exists with that structure and hence the theoretical pre- 
dictions can in principle be compared with experimental 
data concerning the spectrum, susceptibility, etc. 

We study in this paper five different spin ladders char- 
acterized by their local spin S — 1/2, 1 and 3/2 and the 
signs of the exchange of coupling constants along the legs 
Jii and the rungs J±. 
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In the case of the spin 1/2 ladders we discuss the follow- 
ing topics: i) the RVB picture of the antiferromagnetic 
ladder, ii) the equivalence between the ladder state and 
the Haldane state of the spin 1 chain, and iii) the duality 
properties relating the different magnetic structures. 

In the case of the spin 1 ladder we show that the long 
range topological order characteristic of isolated spin 1 
chains disappears and the string correlator decays expo- 
nentially with a finite correlation length. 

The study of the spin 3/2 ladder motivates the def- 
inition of an AKLT state characterized in terms of 9-j 
symbols. 

We compute GS energy densities, spin correlation 
lengths and string order parameters and compare our re- 
sults with those existing in the current literature. 

The organization of the paper is as follows. In section 
II we review the MPM. In section III we particularize the 
MPM to systems which are rotational invariant, where 
the use of group theory leads to a simplification of the 
formalism. In section IV we present our numerical results 
concerning five different spin ladders. In section V we 
summarize our results and present some prospects of our 
work. There are 3 appendices which contain technical 
details or proofs of results presented in the main body of 
the paper. 



II) REVIEW OF THE MATRIX PRODUCT 
METHOD 



Some of the results presented in this section are known 
and can be found in references |]J(|. We also present a 
full account of the formulas and derivations used in refer- 
ence [Q where the MPM was compared with the DMRG 
method in the case of the antiferromagnetic spin 1 chain. 

Let us consider a spin chain or a ladder Bn with open 
boundary conditions, where N denotes the number of 
sites of a chain or the number of rungs of a ladder. To 
describe the low energy properties of -B/v one introduces 
a collection of m states {|a)Ar}™ =1 , which form an or- 
thonormal basis, i.e. N(&\ct')N = Sa.a'- In the DMRG 
these states are the most probable ones to contribute to 
the GS of the superblock -Bjv-i • •-Bjv_i of length 27V, 
formed by adding two sites (or rungs) ••, and a mirror 
image B§_ 1 to the original lattice Bn-i- The basic as- 
sumption of the MPM is that the basis associated with 
B jv and B^-i are related in a simple manner by the 
equation 



\<*)n = J~] A a> p[s] \s) N 

0,s 



)n-i ,N>2 



(1) 



where \s)n (s = 1, . . . , m*) denotes a complete set of m* 
states associated to the N th site (resp. rung) added to 
the chain ( resp. ladder). Eq. (|l]) has to be supplemented 
with the initial data \(3}i- The quantities A a|( g[s] are the 



variational parameters of the MPM, and their determi- 
nation is the central problem one has to solve. This is 
done by the standard variational method. The impor- 
tant point about eq.(Q) is that .A Q)J g[s] does not depend 
on N. Eq.(|l]) is motivated by the truncation method 
used in the DMRG where A a ^[s] dependend on N, i.e. 
A ( ^[s}. When N is large enough one reaches a fixed 



point, i.e. A ( a J[s] 



1 '"' 'i- ->■ Aa.^s], In this manner the ther- 
modynamic limit of the DMRG leads to a translational 
invariant MPM state. 

The condition that both |a)jv and |/3)jv-i form or- 
thonormal basis imposes a normalization condition on 



(2) 



0,8 



It is interesting to count how many variational param- 
eters there are in (|l|). The quantities A Q]| g[s] represent a 
total of m 2 m* variables. We shall assume that all of them 
may be non vanishing. The normalization constraints (|^) 
represent a total of m + m(m — 1)/2 constraints ( m com- 
ing from the diagonal terms a = a' and m{m~ l)/2 com- 
ing from the off-diagonal ones). On the other hand one 
can rotate the basis of states {|a)}™ =1 by an element of 
the orthogonal group 0(m) reducing by m(m — l)/2 the 
number of independent MPM variables. The total num- 
ber of variational degrees of freedom, Na, is then given 
by, 

Na = m 2 rri* — m — 2m(m — l)/2 = m 2 {m* — 1) (3) 

We show in appendix A that the set of A a ,p[s] belongs, 
to the grassmanian manifold, 



0(mm* 



0(m) <g) 0(m(m* — 1)) 



(4) 



As an exercise one can check that the dimension of (^) 
coincides with Na given in (j|). 

In ref. || eq. ([!]) is used to generate an ansatz for the GS 
of periodic chains. In this paper we shall rather use this 
eq. to generate states with open boundary conditions, in 
the spirit of the DMRG. We shall show below that the set 
| a) correspond to ground states with different boundary 
conditions. The use of open boundary conditions leads 
to a simplification of the MPM which is very close to the 
more abstract formalism proposed in |J. 



Correlators of local operators 

Let us use eq.(Q) to compute the expectation values of 
local operators in a recursive way. We shall first consider 
a local operator O n acting at the position n = 1, . . . , TV 
of the lattice. It is easy to get from (|l|) the expectation 
value, 
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\a\O n \a'\ 



N 



T : 



E 



x p \v p ){p p \, (p p \v p/ ) = 5pp> 



(14) 



E/3/3' T aa , M , iv-i^lOnl^OiV-i for n < N 
Y.pO aa , M for n = 7V U 



where 



(6) 



where |i> p ) and are the right and left eigenvectors of 
T with eigenvalue x p {x\ = 1, |t>i) = \v), (pi\ = (p\). 
As a matter of fact all the remaining eigenvalues of T are 
less than one, i.e. \x p \ < 1 Vp 7^ 1. 
In the limit N — > 00 one gets, 



lim Ar ^ oo (a|0i 1 1 ) 0$ • • • OfcV>iV 

<W (p| o( 1 )r ni ~™ 2_1 o( 2 ) ...r n '- i -"'- 1 (9W|w) 



(15) 



aa ', W '=E<^H >WM (s\0\s') (7) 

T can be identified with 1. Eqs.(||,^,0) suggest to in- 
terpret the expectation value isi{o\O n \a')M as a vector 
labeled by the pair aa' , in which case T and O become 
?ti 2 x to 2 matrices. Upon iteration of (J5J) one finds, 

N (a\O n \a') N = T (T N - n d) (8) 

P K /aa',0/3 

More generally, the expectation value of a product of 
local operators is given by, 



N (a\0$0%---0£\a') N 



_ ^ ' frpN — iii ^(l)yii— "2 — 1(^(2) , , jin, — i—n r — 10(r 



(9) 



where AT > n\ > n 2 > ■ • • > n r > 1. The matrix T 
plays a very important role in the MPM. Eqs imply 
that T behaves as a shift operator by one lattice space. 
The basic properties of T follow from the normalization 
condition (^|) which can be expressed as, 



(10) 



which implies that T has a eigenvalue equal to 1. Let us 
call \v) the right eigenvector corresponding to this eigen- 
value. Eq.(|l0|) can be written in matrix notation as, 



The delta function on the r.h.s. of this eq. means that 
the local operators O'"' acting in the bulk, do not modify 
the boundary conditions associated to the various choices 
of a. 

Assuming that T is invertible, one can rewrite eq.([T5|) 
in the following manner, 

lim^oo^HoiV { nl ■ ■ ■ d£ \a') N (16) 
= <W (p\6U(n 1 )6W(n 2 )---0^(n r )\v) 

where 0{n) is defined as 

6(n) = T~ n ~ 1 O T n (17) 

Observe that 1 = 1. The r.h.s. of ( |l6| ) is nothing but a 
spatial ordered product of local operators 0(n), which is 
reminiscent of the radial ordered product that appears in 
•WConformal Field Theory. This connection supports the 
interpretation of T as an euclidean version of the shift 
operator. Under this viewpoint the states \v) and (p\ 
appear as incoming |0) and outgoing vacua (0| that are 
left invariant by the shift operator T. 

We have shown above that the MPM leads in the ther- 
modynamic limit to a sort of discretized field theory char- 
acterized by a shift or spatial transfer operator T and 
local operators (D(n). We can now try to exploit these 
interpretation to extract some physical quantities. 

First of all let us consider the correlator of two opera- 
tors e>W(«i) and 0^(n 2 ). From one has, 



T \v) = \v), v aa , = 5 aa > 



(11) 



( P \ <5« J 1 "! -"2-1 0(2) I ^ 

E^ 1 -" 2 " 1 (p\o^\v p ) ( Pp \d^\v) 



(18) 



On the other hand, let (pi denote the left eigenvector 
of T corresponding to the eigenvalue 1, i.e. 



(P\ T = (Pi <"> ^2p aa 'T aa >.p 0/ = p{3f3> 



(12) 



In the limit when |m — Jia| >> 1 the sum over p 
is dominated by the highest eigenvalue \x p \ of T for 
which the corresponding matrix elements (p\0^ \v p ) and 
{p p \0^ \v) are non zero. If x p < 1 one gets a finite cor- 
relation lcnght £ given by the formula, 



A convenient normalization of (p\ is given by 

(p\v) = 1^J^ Paa = 1 (13) 

a 

For later use we shall diagonalize T as follows, 



i = -\/Hx p \ 



(19) 



In the case where and are both the spin oper- 
ator S, it turns out that the matrix element (p|S|i;) van- 
ishes, and hence the spin-spin correlator is short ranged 
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with a finite spin correlation length £ given by the for- 
mula (|l9|) with \x p \ < 1. The finiteness of £ does indeed 
occur for MP ansatzs which preserve the rotational in- 
variance. However if the latter is broken, as in a Neel like 
state, then £ may become infinite. 

In section III we shall give a formula to compute £ in 
the case of rotational invariant MP ansatzs. 

Another interesting application of ([l6]) is provided by 
the computation of the string order parameter. 

String order parameter 

A spin 1 chain has a long range topological order 
(LRTO) characterized by a non vanishing value of a non 
local operator g(oo) defined as follows JlO|, 

0(00) = lim^oo g(£), (20) 
g(Z) = {S*{t) nti^W S'(0)) 

The AKLT state has 9aklt(oo) = -(2/3) 2 , while 
the spin 1 antiferromagnetic spin chain has g(oo) — 
-0.374325 [0. From eq.© we deduce the following 
expression for (p0|), 

g(£) = {p\T* (e^*)*' 1 &\v) (21) 

In appendix C we show that the operator e™ s " has 
an eigenvalue equal to 1. Denoting by \v st ) and (p st | 
the associated right and left eigenvectors we obtain the 
following expression for the string order parameter 

3(00) = (p\ & \v st ) (p st \ & \v) (22) 

which suggests that g(oo) measures a sort of off-diagonal 
order. 

For antiferromagnetic spin 1 ladders we shall see that 
the LRTO disappears and that the correlator (|o|) is short 
ranged with a finite correlation length £ st . 

Ground state energy density 

Let us suppose we have a translational invariant Hamil- 
tonian of the form, 

N N-l 

H N = "£hU+J2 C+i ( 23 ) 

n— 1 n— 1 

where is an on site (rung) operator while cou- 
ples two nearest neighbour sites (rungs). We define the 
expectation value, 

E% a , — _/v (a\H N \a') N (24) 

which can be computed recursively. From eqs.(^j^) one 
gets 



where h = + h^ 2 \ The hated representation of the 
site hamiltonian ftW is given by eq(R), while the hated 
representation of the hamiltonian hy> is given by, 

= E 77 Vs iV-l,AT<S2Sl|/l^ ) _ l!jv |si4)iV,Ar-l (26) 
X ^a l7 M^ 7l/3 NA^ 7 44]^/H4] 

It should be clear from eqs.(|^, ^6|) which is the hated 
representative of an operator involving an arbitrary num- 
ber of sites. Eq.(p5[) can be conveniently written in ma- 
trix notation as 

\E N )=T\E N ~ l )+h\v), (N>2) (27) 

where \E N ) is a vector with components E^ a ,. Iterating 
(^?j) one gets 

\E N ) = (1 + T + T 2 + ■ ■ ■ + T N - 2 ) h \v) + T™- 1 l^ 1 ) 

(28) 

The geometric series in T can be sumed up and due to 
the eigenvalue equal to 1 it contributes a term propor- 
tional to N, i.e. 

km ±-\E If )=e 0O \v) (29) 

N— >oo IV 

This eq. implies that all the states \o)n have the same 
energy density in the thermodynamic limit, i.e. E^ , = 
5aa' Coo ■ Hence eoo can be identified with the GS energy 
per site for chains or per rung for ladders and it is given 
by, 

Coo = (p\ h\v) =} j p aa ' haa',0/3 (30) 
aa' {3 

This is the quantity one has to minimize respect to the 
MPM parameters. 

The formalism presented above is closely related to the 
DMRG. Eventhough this relation is not the main subject 
of this paper we shall make some remarks ( see H ) . 

MPM versus the DMRG 

Let us suppose that we diagonalize /?, as a m x m ma- 
trix, denoting its eigenvalues as w^, i.e. 

Paa' = W a Saa' (31) 

The eigenvalue eq.(|l^) becomes then, 

£ wlA a p[s] A* a0 ,[s}= 5^w 2 (32) 
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There is a close analogy between eqs(|2|) and (|32|), ex- 
cept for the fact that the order of the labels is exchanged. 
Given the tensor product decomposition s (g) (3 — > a, we 
shall assume that one can reverse the order between the 
states a and f3 in terms of "charge conjugate states" a c 
and /3 C , as follows : s®a c — + (3 C . For example the charge 
conjugate of a state with spin M is another state with 
spin — M. Using this concept we can impose the following 
symmetry condition 



w a A a p[s] = ±wpcApc a c[s], wpc 



wp 



(33) 



which leads to the equivalence between eqs. (||) and (32). 

The relation between the MPM and the DMRG is 
made clear by the construction of the GS of the su- 
per block Bn • Bjy in the following way, 



4>as(3 = Wa A a p[s] 



(34) 



The density matrix that induces ip as p on the block B, 
and which is obtained by tracing over the states in *B^, 
coincides with p aa > — 5 aa > . 



Condition (|33|) guarantees that \ipo) is a state invariant 
under the parity transformation that interchanges the 
blocks B and B^, while leaving invariant the site •. 

It is interesting to observe that the MPM leads to a 
superblock of the form B^ • B^, rather than to the stan- 
dard superblock Bjq • 0. 



Ill) THE MPM APPLIED TO SPIN LADDERS 

In this section we shall apply the MPM to the 2-leg spin 
ladder with a spin S at each site of the chain. The col- 
lection \q)n will be given by the set \JM)n (J < Jmax) 
of states with total spin J and third component M. For 
the sake of simplicity we have only considered one state 
per angular momenta J and M. This will allow us to 
show more clearly the analytic structure of the MPM, 
which can later on be numerically improved by consider- 
ing multiplicity. This has already been done in the case 
of spin chains in references 

The states added at each step of the MPM are the ones 
that appear in the tensor product decomposition of two 
spin S irreps, i.e. S (g> S = 0®l(B ■ ■ ■ ®2S. These states 
are labelled by |A/x) where A = 0, • • • , 2S is the total spin 
and fi = —A, . . . , A is its third component. 

Using these notations we propose the following recur- 
rence relation for the states \JM)m% 



where 



\JxM l ) N = Y J A) l:h |(AJ 2 ), J\Mx) 

XJ 2 



|(AJ 2 ), JiMOat 



N 



(35) 



(36) 



(V, J 2 Af 2 |AJ 2 , JiMx) \\n) N ® \JiM- 



2 N-l 



In (f36J) the quantity (Xfj,, J 2 M 2 |AJ 2 , JiM x ) is the 
Clebsch-Gordan coefficient corresponding to the decom- 
position A <8> J 2 — * J\- Comparing eqs.(Q) and ( |35| ) 
we obtain the following relation between the symbols 
A,j 1 m 1 j 2 Af 2 [A£t] and the rotational invariant symbols 

Aj iMi .j 2 m 2 [M = A x hJ2 (A/x, J 2 M 2 \ AJ 2 , JxMx) (37) 

The use of rotational invariant basis reduces consider- 
ably the number of independent variational parameters 
and consequently increases the power of the MPM 1,0. 

The variational parameters Ajj are subject to the 
CG condition, 



A JiJ 2 



unless | A - J 2 \ < Ji < | A + J 2 | 



(38) 



Using (|37|) and the orthogonality properties of the CG 
coefficients, the normalization conditions (0) become, 



A,,7 2 



JiJ 2 \ 



1, VJi 



(39) 



At this point we can just take eq.(|37|) and plug it into 
the corresponding formulas of section II in order to derive 
expectation values, the GS energy density, etc in terms 
of A j j 2 . There is however a more efficient way to do this 
by using group theory. The application of the Wigner- 
Eckart theorem will allow us to express all the results 
in terms of reduced matrix elements of the operators in- 
volved as well as the 6j-symbols. In our derivations we 
shall follow the same steps as in section II, leaving the 
technical details to appendix B. 



Correlators of Invariant Tensors 

Let us denote by 0^ k ' an irreducible tensor of total 
angular momentum k, whose components are labeled by 

(k) 

O m ,M = —k, . . . , k. The spin operators S correspond 
to k = 1. Let us suppose we have two irreducible tensors 
with the same total angular momenta k, 0^ k ' A ^{n) and 
CK fc ' S ) (rn), acting at the positions n and m ( N > n > 
to > 1) of the ladder. The scalar product of these two 
operators is defined as 



(40) 



-M n (k,A) 



(n).O^(m) 



The basic result we derive in appendix B is, 
N (JxM\ 0^ A \n) ■ 0^ B \m) \JxM) N 



(41) 
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ErpN-n Q{k,A) /jro— m-l\ Q(k,B) 
J1J2 J%,J&J4 V k I J 3 J i} J 5 J 6 J$J6,Jy 



where T and T^> are defined as, 

T Ji,J 2 = J2\ (A/iJ 2 ) A \j 2 

Jl J2, J3 J A 

x (-l) A+fc+Jl+J V(2^i + 1)(2J 2 + 1) {jIjII 



Using (B5|) we find 



■h - E ( T *.* + *JlJ, ) . > 2) (- r .i)) 

J 2 



(42) 



(43) 



while and O^ 5 ) are defined in Appendix B. 

Eq.(^) is the invariant version of (^J) involving only 
two operators. In order to obtain the thermodynamic 
properties of ( [4l| ) we use the properties of the transfer 
operator T. The normalization conditions ( |39| ) imply the 
following conditions on T, 



J 2 



(44) 



Let us call pj the left eigenvector of Tj lt j 2 with eigen- 
value 1, i.e. 



(45) 



Ji 



Using eqs.([l4, 45) into (^Tj) and taking N >> 1 we get 
lim^oo iv<JiM| G( fe ^)(n) • 0( fc ' S )(m) | JxM)^ (46) 

pn 
l k 



where we use a matrix notation in J-space with the con- 
vention Vj — 1, VJ. From eq.(^) we deduce that the 
correlation length associated to the scalar product of two 
irreducible operators with angular momentum k, is given 
by the highest eigenvalue of the matrix defined in (43). 
The spin-spin correlation length is obtained by looking 
at the highest absolute eigenvalue of T\. 



Ground state energy density 

The Hamiltonian of the 2-leg ladder has the form pro- 
posed in ( p3|) where is the rung Hamiltonian and 
is the leg Hamiltonian, 



hn = J± Si(n) • S 2 (n) 



(47) 



^S+i^ll (S 1 (n)-S 1 (n + 1) + S 2 (n)-S 2 (n+1)) (48) 

S a (n) is a spin S operator acting on the n = 1, • • • , N 
rung and the a = 1, 2 leg of the ladder. 

As in ([m]) we define the expectation value of the ladder 
Hamiltonian, 



E'J = N (JM\H N \JM) 



N 



(49) 



where h = + (/i« and can be found in 
appendix B). 

Iterating eq. ( |50|) and using the properties of the matrix 
T we can inmediately get the large N limit of the energy 
(1), 



1 



lim —E 1 ? = e™, V J 



N->oo N 



J c oo ? 



where the GS energy density is given by, 

eoo = (p\ h \v) = PJi h JlJ2 
J1J2 



(51) 



(52) 



At this point let us summarize the main steps of the 
MP algorithm hereby proposed, 



• Solve the normalization conditions (|39j) expressing 
A j j in terms of a set of linearly independent vari- 
ational parameters. 

• Find the eigenvector pj of the matrix T. 

• Minimize the GS energy density (|52|) with respect 
to the independent variational parameters. 

We will now comment on how these three steps can be 
implemented. 



Solution of the normalization conditions 



We shall suppose in the rest of the paper that the pa- 
rameters A j j 2 are all real. Hence the normalization con- 
ditions 



{A) iJ2 f = l, V Ji 



(53) 



imply that the set {Aj j } for J\ fixed are the coordinates 
of a sphere whose dimension depends on the allowed val- 
ues of J and the CG conditions (|38|). Let us call A™ 3 * 
the highest coordinate, in absolute value, i.e. 

^r = ^L suchthat I^XI^I^I VA,J 2 (54) 

If AJ™ > 0( resp. Af^ < 0) we can think of it as the 
north ( resp. south) pole of a sphere, whose neaghbour- 
hood can be described by the stereographic coordinates, 



J JlJ2 



= A 



JlJ: 



IK 



■J1J2 



< 1 



(55) 



Notice that Xj° Lq = 1. The remaining coordinates are 
the independent variational parameters used in the mini- 
mization of the GS energy. The solution of the constraint 
(p>3|) finally reads, 
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-1/2 



A JlJ2 - X JlJ2 ^2 \ X JlJ 2 



, ej = ±l 



,A' J' 



(56) 



where ej t = 1(— 1) corresponding to the north ( south ) 
polo of the above mentioned sphere. 



Determination of p.j 

The solution of the eigenvalue problem of eq.([l5"|) can 
be done numerically. However for a ladder with spin 
S = 1/2 it can also be solved analytically which will 
allow us to make some considerations on the nature of 
pj. In the case where 5=1/2 the allowed values for A 
are and 1. Hence the unique non- vanishing entries of 

A Jijfa &Te A °JJ> A JJ> A JJ+i and A JJ-i- Similarly from 
eq.(|42|) the non- vanishing entries of T are Tj t j,Tj j+i 
and Tj ; j-i. The set of eqs. we have therefore to solve 
read explicitely as, 



Tj : j + Tj : j + i + Tjj^i — 1 
Pj Tjj + pj + i Tj + i,,/ + pj_iTj_i t j = pj 

Ej=o x Pj = 1 

The solution of these eqs. is given by, 

uj 



(57) 



Pj 



where 



«:: I- "J \[ ~T~~~~~~ (J ■-'») 



L=0 



(58) 



(59) 



we are assuming that J = 0, • • • , J m ax- 

Eqs.(^8|, [59| ) imply that pj is always positive, in agree- 
ment with the Perron-Frobenius theorem applied to the 
matrix T, whose entries are all non-negative. In ref. JtJ it 
was shown that the values of pj are intimately related to 
the eigenvalues of the density matrix that appear in the 
DMRG. These and other facts suggest that the MPM 
is in fact equivalent to the DMRG, specially when the 
number of states kept m becomes large. 

This completes the presentation of the formalism. 



IV) NUMERICAL RESULTS 

In this section we shall apply the MPM to five differ- 
ent spin ladders, corresponding to different choices of the 
spin 5 and signs of the coupling constants Jm and J±. 
We shall denote every of these ladders as AA$, AF$, and 
FAs where A and F stands for antiferromagnetic or fer- 
romagnetic couplings. Thus for example AFs is a spin 



S ladder with antiferromagnetic couplings along the legs 
and ferromagnetic couplings along the rungs. With these 
notations we will study below the following cases: AAi/ 2 , 
AFx/2, FA1/2, AA\ and AA^/2- Within each case we will 
highlight a particular aspect, which the MPM helps to 
clarify. 



yL4 1/2 -ladder: The dimer-RVB state 

This is the most studied spin ladder. Its properties are 
well known and can be summarized as follows. In the 
weak coupling regime, i.e. Jj_ << Jm, the gapless spin 
1/2 chains become massive by the interchain coupling 
which is a relevant operator of dimension 1 |l3Hl5|] . The 
magnitude of the gap is proportional to Jj_ . In the inter- 
mediate coupling regime, i.e. Jj_ ~ J11 the spin ladder can 
be mapped into the 0(3) non linear sigma model (NLSM) 
with no topological term [p~6| — P~S|| . This model is known to 
have a spin gap. From numerical studies the magnitude 
of the spin gap A and the spin correlation length, in the 
isotropic case J± = J11 = J, are given by A = 0.502J 
and £ = 3.2 respectively ]l9|"p2]|. In the strong coupling 
regime J± >> Jm, the most appropiate physical picture 
of the GS and excitations is given by the RVB scenario 
proposed in pc| , and supported by DMRG p0[ , mean 
field |^3| and variational calculations |24|]. In the latter 
work a recurrent variational ansatz (RVA) was proposed 
to generate the dimer-RVB and generalizations of it. The 
RVA method is a MPM based on 2 nd and higher order 
recurrent relations, while the standard MPM is based on 
a 1 st order relation. We shall see below that the MPM 
applied to ladders essentially contains the RVA, and that 
the numerical results are improved. 

Let us first consider the case where the MPM states 
\JM)m are choosen to be a singlet and a triplet, i.e.. 
J = and 1. In this case eq. ( p5|) is depicted in fig. 1. 
There are a priori 5 non vanishing MP parameters sub- 
jected to 2 normalization constraints, leaving a total of 3 
independent parameters. 

Fig. 2 shows A JiJ2 as functions of x = J\\/J±. In the 
whole range of coupling constants the most important 
amplitudes are A^ Q and A\ while is essentially zero. 
The latter amplitude correspond to having only a single 
bond among two rungs ( see fig. lL which is forbidden 
in the dimer-RVB picture of refs. (2(J and [p4| . 

In table 1 we show the GS energy density obtained 
with the MPM for J max = 1 and 2, together with the 
RVA, mean field and Lanczos results. In table 2 we give 
the spin correlation length computed with the MPM and 
the RVA. 

There is an appreciable improvement in the numerical 
results of the MPM respect to the RVA, specially for the 
spin correlation length. 
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J\\/Jl. ^raax — 1 


Jmax — 2 RVA 


Mean Field Lanczos 


0.0 


375000 

VJ . kj I <J V / VJ U 


0.375000 0.375000 


375000 




0.2 


383199 


0.383199 0.383195 


382548 




0.4 


40Q607 


0.409608 0.409442 


405430 

VJ . t:UU t:OVJ 




0.6 


453509 


0.453513 0.45252 


442424 




0.8 


0.510504 


0.510523 0.507909 


0.489552 




1.0 


0.575924 


0.575970 0.571314 


0.542848 


0.578 


1.25 


0.664776 


0.664867 0.657551 


0.614473 


0.6687 


1.66 


0.819656 


0.819834 0.808438 


0.738360 


0.8333 


2.5 


1.152056 


1.152416 1.13384 


1.002856 


1.18 


5.0 


2.172002 


2.172878 2.13608 




2.265 



Table 1. GS energy per site —e 00 /2J± of the ladder 
AAi/2- The first two columns are the MPM results. 

The RVA results are obtained with a third order 
recursion formula (24|]. The mean field and Lanczos 
results have been obtained in references |^3| and p5[ | 
respectively. 



J\\/ J 1. <^max 1 


T — 2 

'-'max ^ 


RVA 


0.0 


0.00000 


0.0000 


0.000000 


0.2 


0.5300 


0.5303 


0.437166 


0.4 


0.8057 


0.8081 


0.608323 


0.6 


1.0652 


1.0740 


0.751286 


0.8 


1.2753 


1.2945 


0.866958 


1.0 


1.4282 


1.4593 


0.959249 


1.25 


1.5572 


1.6018 


1.04877 


1.66 


1.6802 


1.7413 


1.15205 


2.5 


1.7903 


1.8698 


1.26951 


5.0 


1.8747 


1.9711 


1.38532 



Table 2 Spin correlation length of the ladder AA 1 / 2 - 
The first two columns are the MPM results. The RVA 
results are those of ref |24l . 



AFi/2-ladder: Relation with the spin 1 chain 

The ladder with magnetic structure AF is interesting 
because it is intimately related to the spin 1 chain p6| . 
This relation can be clearly seen in the strong coupling 
limit — J± >> Jii, since it leads to an effective spin 1 on 
every rung, which are effectively coupled antiferromag- 
netically along the legs. The effective Hamiltonian can 
be derived from (|7|,[48|) and reads Q , 



Tj ladder 



J±]N- 



1 



Jii 



E 



Scff • s cff 



(60) 



where S c ff(n) = Si(n) + S-2(n) is the spin 1 operator 
acting on the n th rung. The term proportional to J± 
comes from the rung Hamiltonian when diagonalized in 
the spin 1 sector. This eq. implies the following relation 
between the energies per site of the AF-ladder and the 
spin 1 chain, 



1 , , 1 

-8^ + 4 J '' 



„cff 



(61) 



In table 3 we give the GS energies of the ladder 
parametrized in terms of the effective energy e^. We 
also give the spin correlation length. We have made two 
choices of MP states | JM ) m . One for which J = and 1 
and the other for which J — 1/2 and 3/2. The latter one 
corresponds to having a single spin 1/2 at the boundary 
of the ladder. These two choices have an analogue for 
the spin 1 chain. For integer J's the MP parameters do 
not vary in the whole interval < Jn < 1.66 and take 
the values, 



A° 



A 1 

^10 



A° 1X = 0.000,4^ = 1.000 
= -0.577, ^ = 0.816 



(62) 



For half- integer J's the MP parameters do not vary in 
the whole interval < Jii < 5 and take the values, 



Al 1 = A% 3 = 0.000, A\ 



0.989 



(63) 



A\, =0.148, ,4! i 



-0.953, A! 3 



-0.303 



Note that for half integers J there is one more vari- 
ational parameter. The most important amplitudes in- 
deed correspond to the formation of an AKLT state with 
a single bond connecting every effective spin 1. 



-J\\/J± 



-e 



ctt ,int 



TiiEEaTT 



FKaTT 



(0.0, 1.66) 1.333333 1.399659 0.910 2.5997 



2.5 



1.363970 1.399659 1.9682 2.5997 



5.0 



1.498539 1.399659 1.9607 2.5997 



Table 3. The exact values of the GS energy density of a 
spin 1 chain and its correlation length are given by 
eoo = -1.4014845 and £ = 6.03 (n). 

The values of and £ in the integer J case coin- 
cide with those of the AKLT state, while the ones of 
the half-integer J case coincide with those obtained with 
a MPM method applied to the spin 1 chain Q| where 
one keeps two MPM states with J = 1/2 and 3/2. These 
results provide additional support for the equivalence be- 
tween the AFi/2 ladder and the spin 1 chain in the strong 
and intermediate coupling regimes observed previously 
by other authors [f26fl. 



FA 1/2 -ladder 

In the strong coupling regime the ladders FA 1 / 2 and 
AF1/2 have similar GS energies and correlation lengths 
( see tables 1, 2 and 4). The MP parameters display 
also a similar behaviour although some of them are in- 
terchanged ( see figures 2 and 3). The physical reason of 
this is the common GS in the case where J\\ = 0, given by 
the coherent superposition of valence bonds in the rungs. 
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-J\\/J± 


-eoo/2J_L 




0.0 


0.375000 


0.0000 


0.2 


0.381754 


0.5140 


0.4 


0.399295 


0.7577 


0.6 


0.424396 


1.010 


0.8 


0.454891 


1.277 


1.0 


0.489324 


1.554 


1.25 


0.536374 


1.895 


1.66 


0.619895 


2.381 


2.5 


0.803434 


2.992 


5.0 


1.376973 


3.520 



Table 4. GS energy per site and correlation length of 
the ladder FA 1 / 2 . 

The relation between FA 1 / 2 and AF 1 / 2 is part of a 
more general relation involving also the ladder AA 1 / 2 and 
can be established by means of a type of transformations 
called dualities in ref. If29|]. 



Duality properties of spin ladders 

On a 2-leg ladder one can define 3 types of dualities 
called U, T and S, which mix or leave invariant the lad- 
der's magnetic structures AA, AF and FA p9| . 

The U duality maps a Hamiltonian with couplings 
constants Jy , J_l into a ladder with couplings constants 
JF , JV where, 
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: J,, (S 1 (n)-S 1 (n + l))/(S 1 (n)-S 2 (n 



1)) (64) 



Under U the leg-bonds are transformed into diagonal 
ones while the rung-bonds are left invariant. The signs of 
(Si(n) • Si(ra+ 1)) and (Si(n) • S 2 (n + 1)) are determined 
by those of Jn and J± respectively. Thus U acts on the 
magnetic structures as follows, 



AA 
J,, (AA) > 
J AAA) > 

AF 
J\\(AF)>0 
J±{AF) < 



FA 

JF(FA) < 
Jf(FA) > 



AF 
JF(AF) > 

jI(af) < o 



(65) 



(66) 



In fig. 4 we show JF (FA) and JF (AF ) as functions of 
Jn (AA) and Jm (AF) respectively. 

The GS energy density of the ladder with coupling con- 
stants JF, Jf i s a lower bound of the original GS energy 
il, i.e. 



3 (J,f (FA), J] 7 (FA)) < eoo (J||(AA), J±(AA)) 
3 (J^(AF), Jf (AF)) < eoo (J||(AF), J ± (AF)) 



(67) 



In fig. 5 we show the validity of these inequalities, 
which in the strong coupling limit almost become identi- 
ties. In fig. 6 we show the correlation legths for both AA 
and the transformed FA ladders. Again in the strong 
coupling limit they become very close. 

The T transformation consists in the replacement of 
the vertical bonds by diagonal ones, i.e. 



J( = J,, 

Jl = J ± (Si(ra) • Si(n + l))/(S 1 (n) • S 2 (n + 1)) 



(68) 



which leads to the following action on magnetic struc- 
tures, 



AA 
J|| (AA) > 
Ji(AA) > 

FA 
J\\(FA)<0 
J± (FA) > 



AF 



Jf(AF) > 
Jl(AF) < 



FA 
JT(FA) < 
J±(FA) > 



(69) 



(70) 



In fig. 7 we plot the energies associated to the FA lad- 
der and its T transformed, which satisfies the inequality, 

e oc (Jf(FA),Jl(FA)) < e oc (J ]l (FA),J±(FA)) (71) 

The convergence of both curves in the weak coupling is 
in agreement with the bosonization arguments employed 
in §9). 

Finally the S transformation is defined by the replace- 
ment of vertical bonds by horizontal ones and viceversa, 



J,f = I J x (Si(n) ■ S a (»))/(Si(n) • S 1 (n 



Jf = 2 J,, (Si(n) -Si(n- 



1)> 

l))/<S!(n).S 2 (n)) 



(72) 



The factors 2 and 1 /2 are explained by the fact that there 
are two leg-bonds for each rung-bond. Eqs.(^2h imply, 



AF 
J,, (AF) > 
Jl(AF) < 

AA 
J,l (AA) > 
Jx(AA) > 



FA 
J[(FA) < 
Jf (FA) > 

AA 
jf(AA) > 
Jf (AA) > 



(73) 



(74) 



In fig. 8 we plot the energies of the AA ladder and its 
transformed which satisfy the inequality, 

e oc (J | f (AA), Jf (AA)) < eoo (J||(AA), J ± (AA)) (75) 

Note that in the region Jn ~ J± both energies get very 
closed. Fig. 9 shows the spin correlation length for the 
AA ladder and its S transformed, displaying the same 
pattern as fig. 8. In summary we have found further 
numerical evidence of the duality properties of the 2-leg 
ladder proposed in p9[). 
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j4j4i-ladder: Short-range string order 

In table 5 we show the GS energy density and the spin- 
correlation length of the ladder AA\. Observe that the 
correlation length is longer that the one of the spin 1/2 
ladder. 



J\\/Jx 


-eoo/2Jj_ 


£ 


0.0 


1.000000 


.00000 


0.2 


1.055719 


1.0114 


0.4 


1.206557 


1.8318 


0.6 


1.407358 


2.3852 


0.8 


1.631166 


2.6762 


1.0 


1.867327 


2.8227 


1.25 


2.172905 


2.9042 


1.66 


2.688880 


2.9286 



Table 5. GS energy per site and correlation length of 
the AAi ladder. 

As mentioned in the introduction a spin 1 chain has a 
long range topological order (LRTO) characterized by a 
non- vanishing 5(00). In appendix C we give an analytical 
expression for g(oo) in terms of the MP parameters of the 
spin 1 chain. 

However when two spin 1 chains are coupled antiferro- 
magnetically the LRTO disappears and the string order 
parameter g(l) decays exponentially as e~ i '^ . We call 
£ st the string correlation length, and its value together 
with the spin correlation length are shown in fig. 10 as 
functions of the ratio J\\/J±. In the weak coupling limit 
where J11 / Jj_ — * 00 we expect £ st to diverge, recovering in 
that way the LRTO of the uncoupled chains. The value 
of £ st is obtained by the formula (|l^) with x p the highest 

eigenvalue of the operator e l7rS i (see appendix C). 

An intuitive way to understand the breaking of the 
LRTO is given by the AKLT picture of ref @. An AKLT 
state is a valence bond state where every spin 1 is repre- 
sented as a symmetrized product of two spins 1/2, and 
such that every of these "elementary" spins is linked by 
a bond to one of the spins 1/2 on its neighbours. In 
this way all the spins of the chain are connected by a 
sucession of nearest neighbour links. When we couple 
antiferromagnetically two spin 1 chains there is the pos- 
sibility that two parallel bonds along the legs become 
two parallel bonds along the rungs as shown in fig. 11. 
Thus the two infinite parallel arrays of connected bonds, 
characteristic of the uncoupled chains, effectively breaks 
into a collection of fluctuating islands whose size is of the 
order of £ st . Everyone of these islands is a sort of closed 
spin 1 chain ( fig. 11). 

The finite value of £ st at the origin of fig. 10 is due to the 
fact that e l7rSl has indeed a finite value when computed 
on the singlet formed by two spins 1 on a rung, 

(e^) rU n g = J2 (-l) m (lml-m|00) 2 = -i (76) 

m=±l,0 '* 



which leads to £ st (J|| — 0) — l/ln3. Fig. 10 suggests the 
existence of 3 different regimes. In the weak coupling 
regime where £ st > £ the ladder can be effectively con- 
sidered as a collection of weakly interacting closed spin 
1 chains. In the strong coupling regime where £ st > £ 
the bonds are meanly along the rungs and the interbond 
coupling is small. Finally there is an intermediate region, 
with £ st < £, where the islands of spins interact strongly 
with their neighbours. 

The AA 3/2 -ladder 

In table 6 we give the GS energy densities and spin cor- 
relation lengths of the ladder AA 3 / 2 - As one may expect 
the correlation length is longer than for the spin 1 and 
1/2 ladders. This fact agrees with the results obtained 
by mapping the spin ladders into the NLSM |To| Q . 



J\\/J± 


-e ao /2J± 




0.0 


1.875000 


0.0000 


0.2 


2.054760 


1.8760 


0.4 


2.449827 


3.3099 


0.6 


2.911353 


3.9475 


0.8 


3.400562 


4.2401 


1.0 


3.904988 


4.3829 


1.25 


4.548607 


4.4624 


1.66 


5.623131 


4.4900 



Table 6. GS energy per site and correlation length of 
the AA 3 / 2 ladder. 

AKLT states for ladders 

The spin 3/2 2-ladder offers the possibility of con- 
structing an AKLT state with a valence bond connect- 
ing every spin 3/2 to its three nearest neighbours. More 
generally, let us consider a ladder with spin 5 > 3/2 and 
three integers p,q,r > 1 satisfying the eq. 25 = p + q + r. 
Then one can define an AKLT state, denoted by the 
triplet (p, q, r), by linking the 25 "elementary spinors" 
of each spin to the ones in its neighbours following the 
pattern shown in fig. (12). The AKLT states (p,q,r) and 
(q, p, r) when p ^ q correspond to dimerized ladders and 
they differ by the translation of one unit space along the 
legs. 

The spin 3/2 AKLT ladder corresponds in the above 
notation to (1,1,1). This state contains in fact a spin 
and a spin 1 state which can be generated by the MP 
equation (|35|) where the amplitudes Ajj are given by 
9-j symbols, 

f 1/2 1/2 J 2 ) 

A x JiJ2 = 3 V(2J 2 + 1)(2A + 1) \ 3/2 3/2 A \ (77) 
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In this eq. J\, J2 = and 1, while A = 0, 1, 2. 

The proof of (|77]) follows from the definition of the 9-j 
symbols as the coefficients that give the change of basis 
when coupling in two different ways 4 angular momenta, 
namely Q 

^{3ih{Jxz)hh{J2i)J) (78) 
= J2 V(2Ji2 + 1)(2J 34 + l)(2Ji 3 + 1)(2J 24 + 1) 

!3x h J12 1 
33 k J34 > i>{3i3i{Ji2)h3i{Jai)J) 
J\3 J24 J ) 

where ip {3133^13)3234^24) J) is a state with angular mo- 
mentum J obtained by the tensor product decomposition 
J13 <%> J24 — > J, which in turn are obtained by the decom- 
positions ji ® j 3 — > J13 and j 2 ® 34 — > J24- 

One may check that the normalization conditions j39| ) 
holds for ( |77j ) , as a consequence of the orthogonality con- 
ditions satisfied by the 9-j symbols j30|. The GS energy 
per site and the spin correlation length of the AKLT state 
( f77| ) in the case where J\\ = Jj_ = J are given by, 

e AKury 2J = _ 3 263536, £ AKLT = 1.116221 (79) 

This state has a much shorter correlation length than 
the MP state that minimizes the GS energy of the AA 3 / 2 ~ 
ladder (see table 6). The GS energies of both states are 
also quite different. We conclude from these facts that 
the spin 3/2 AKLT state does not give a good description 
of the GS of the AA 3/2 -ladder. 

A generic AKLT state of the type (p, q, r) when p 7^ q 
has to be described by alternating MP amplitudes de- 
pending on the eveness of the site. Thus for even sites 
one has 

r § § j 2 

=(q + r + l) V(2 J 2 + 1)(2A + 1) I S S X 

q+r q+r j 
(.2 2 Jl 

(80) 

where J\ = 0, . . . , q; J% — 0, . . . , p and A = 0, . . . , 2S — 
r. For odd sites the corresponding MP amplitudes are 
obtained by interchanging p and q in (|80|) . 

V) CONCLUSIONS AND PROSPECTS 

Let us summarize the main results obtained in this 
paper. 

• We have presented a rotational invariant formula- 
tion of the MPM which allow us to express the GS 
energy density, the correlation length and the string 
order parameter, in terms of invariant objects. This 
reduces considerably the number of independent 
MP parameters used in the minimization process. 



• We have improved the numerical results concerning 
the GS energy density and spin correlation length 
obtained previously with other approximate meth- 
ods as those of references |^J,^3|. The considera- 
tion of MP ansatzs with multiple states per spin 
will certainly lead to better results. 

• We have shown the equivalence between the ladder 
AF 1 / 2 and the spin 1 antiferromagnetic Heisenberg 
chain. The MPM applied to both systems shows 
strong numerical coincidences for the GS energy 
and correlation length. This agrees with the results 
obtained previously by other methods |2(] |2g] 

• We have found numerical evidences for the duality 
properties proposed in |^9| for the spin ladders with 
magnetic structures AA, AF and FA. 

• We have shown that there is a breaking of the long 
range topological order of the spin 1 chains when 
they are coupled in a 2 legged ladder. A physical 
picture of the GS of the spin 1 ladder is given in 
terms of resonating closed spin 1 chains. 

• We have constructed AKLT states for 2 legged lad- 
ders with spin S > 3/2, showing that the corre- 
sponding MP parameters are given by 9-j symbols. 

• We have suggested a relation between the MPM 
and the DMRG based on the density matrix that 
appear in both methods (see also 0). We conjec- 
ture that the minimization of the GS energy eoo 
can be transformed into an eigenvalue problem on 
a super block Bn • B§. 

In summary we have shown the adequacy of the MPM 
to study the 2 legged ladder, specially in the strong and 
intermediate coupling regimes. This is made possible 
from the fact that these ladders are finitely correlated. 
Hence one may expect that even spin ladders with a finite 
number of legs could be described by the same technique, 
although with a larger number of states m. On the other 
hand odd legged ladders are not finitely correlated and 
they cannot be properly described in the large N limit 
within the actual formulation of the MPM. An interesting 
problem is the application of the MPM to 2D systems, 
which can be thougth of as ladders with a large number 
of legs. It is clear that one should choose a collection of 
the most representative states for the rungs to be added 
after each iteration of the MP recurrence equation. 
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APPENDIX A: THE MP ANSATZ AND THE 
GRASSMANNIAN MANIFOLDS 

In this appendix we shall give a proof of eq. (||) which 
gives a precise mathematical meaning of the coefficients 
Aq,/3[s] defining a generic MP ansatz. 

In the r.h.s. of eq. (Q) we have a generic vector of dimen- 
sion n = mm* while on its l.h.s. the vector has dimension 
m. Hence eq.Q amounts to a choice of a TO-dimensional 
linear subspace of R n in the case of A a p [s] real or a com- 
plex subspace of C™ in the case of j4 Q( 3[s] complex. Let 
us call the set of all these subspaces as Af„, m (R) and 
M„ !m (C) for A Q/ 3[s] real and complex respectively. The 
group 0(n) ( resp. U(n) ) acts transitively on M„ jTO (R) 
( resp. M n ^ m (C)), which leads to the result |3l"|] 



M„, m (R) 



0(n) 



0(m) ® 0(n — m) 



M n>m (C) = 



U(n) 



U(m) (g) U(n — m) 



(81) 



(82) 



In ( pl|) the groups 0(m) and 0(n — m) are identified 
with the subgroups of 0{n) consisting of those elements 
leaving fixed every vector of a given (n — m)-dimensional 
subspace and of its orthogonal complement, respectively. 
Similar arguments lead to eq. (|S2|) . Af nim (R)(M njm (C)) 
are called the real (complex) Grassmannian manifolds. 
Taking n = mm* in (|^) we get eq.(^). 

As a simple illustration of these eqs. let us consider the 
case of a MP ansatz that generates a single state \GS) n 
(to = 1), i.e. 

\GS)s = ^2A[8]\s)s9\GS) N -i (83) 

S 

with A[s] G R. The normalization condition ^ reads, 



(84) 



Thus A[s] belongs to the (m* — l)-dimensional sphere 
SO(m*)/SO{m* - 1). Upon the identification of A[s] 
and — A[s] we get the (m* — l)-real proyective space 
M TO »,i(R) = SO{m*)/SO{m* - 1) ® Z 2 . 
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APPENDIX B: THE ROTATIONAL INVARIANT 
MPM 

Group theoretical preliminaries 

Before we give the proof of the main formulas of section 
III we shall review some basic definitions and results in 
group theory Q . 

An irreducible tensor with angular momentum k is an 
operator T£p (M — k,---,—k) which satisfies the fol- 
lowing commutation relations with the total angular mo- 
mentum operator J, 



[J Z) T$]=MT$ 

[J x ± % Jy,T$] = y/k(k+l)-M(M±l) T$ ±1 



(85) 



The scalar product of two irreducible tensors and 
U< fc ) with the same spin k is defined by, 



^ (-l)- M T$U«l (86) 



M=-fc 



The Wigner-Eckart theorem reads, 
(JM\T^ k) \J'M') 

= (-iy- M (_ J M I J M .) (jiit^iuo 



(87) 



where the 3-j symbol is related to the CG coefficient by 



J k J' 

-M fj, M' 



(-1) 



J-k-M' 



V2J> 



(J - Mkfj,\J' - M') 



(88) 



The quantity ( J\ \ | J') in © is called the reduced 
matrix element of the operator TW . As an example we 
give the reduced matrix element of the spin operator S, 



(S\\S\\S) = y/S(S+l)(2S + l) 



(89) 



Let \0c1j1a2j2JM} be a state with total angular mo- 
menta J and third component M, appearing in the ten- 
sor product decomposition (aiji) <8> (0^232), where (aj) 
denotes a state with total angular momentum j and a 
labels other possible quantum numbers. We shall need 
below the following results. 



(aijia 2 j2JM\(T [ { 



a'JMJ'M') 



x(a m ||T«||aiJi) (a 2 j 2 \\T^\\a' 2 f 2 ) 



{aijia 2 hJ\\TW\\a[j[a 2 j 2 J>) 
<) ,<)■ •/ (-\\ii+h+J'+k J 3\ J h 



(90) 



(91) 



x V /(2J+1)(2J' + 1) (a^WT^Wa'^) 



x vW + W+1) (a2j2||T 2 fc) ||a 2 j 2 ) 



(92) 



The subindices 1 and 2 in and T 2 K> mean that 
the corresponding operators acts on the states labelled 
as (aiji) and («2j2) respectively. 

Recursion relations for the scalar product of in- 
variant tensors 

We want to prove eq.([i"l"|). 

Using eq.(pq) we easily get for N > n > m, 



-.(*) 



n(JiM\ {k ' A) {n) ■ {k ' B \m) \.hM) N (93) 

= E T -ruJ, N-l{JiM\ 0^ A \n)-0^ B \m) \J 2 M) N - X 
J2 



where Tj ly j 2 is given in fl42|). Iterating (|93|) N — n times 
we reach the situation where N = n. This produces the 
term T T ~j n in @. Next we need to compute the matrix 
element, 



,{JiM\ (k ' A \n) ■ {k > B \m) \J X M) 7 



(94) 



E 



A) 2 T ) A) 3 , 

J1J2 / J1J3 



x„((A 2 J 2 ), JiM\ 0^ A \n) ■ 0^ B \m) |(A 3 J 3 ), JiM)„ 

The matrix element on the r.h.s. of (0) has the form 
described in (^p|), which yields, 



,((A a J a ), JxM| 0^)(n) • G^ s )(m) |(A 3 J 3 ), JiM)„ (95) 

- r_iv / i+^+A 3 / A2 J 2 J l 
1 ' 1 J 3 A 3 fc 

x „(A2||0( fe ^H||A 3 )„ „_ 1 (J 2 ||0( fe . B )(m)||J 3 )„- 1 
Introducing ( |95| ) into ( p4] ) we find 

„(J!M| 0( fe - A >(n) • 0^ s )(m) \J x M) n (96) 

= E n-i(J a ||0 (fc ' B) M||J 3 )^i 
J2J3 



where 



Ojjj = Ea a ( A j 2 j A ) 3 



(97) 



x(-l)^+* £ ^}(A 2 ||0(^)||A 3 ) 



The next step is to apply the MP ansatz ( p5| ) to 

„(Ji||0( fc ^(m)||J 2 )„ (98) 



13 



E * 



A 2 

JiJa 



where 



X„((AlJ 3 ), Jl|| (M3) M IKA2J4), J 2 )r 

For n > m we can use (p^), getting 

n ((AiJ 3 ),Ji||0 (fc,B) (m) ||(A 2 J 4 ),J 2 )„ 



(99) 



V 1>Ja - ^ E Q A ( A + !) - S ( 5 + !)) I A U I' 

(105) 

Jl,.74 II / / V Ji J 2 J2J4/ J1J3 J3 

J4,Ai ,...,A4 



<$AiA 2 (-1) 



Ai+.7 1 + J 4 +fe 



V(2Ji + 1)(2J 2 + 1) 



v / 1 \ I + A3+A4 fA 2 Ai fA 3 A 4 



(106) 



and 



= (-l) Jl+J3 v / (2J 1 + l)(2A 1 + l)(2A 2 + l) 



Plugging (99) into ( Pq ) we get, 

n (J 1 ||0( fc ' i3 )(m)||J 2 )„ 



(100) 



Xv /5(5 + l)(25 + l) 





A 2 






A 2 






J2 






S 


1} 



E ( T k)jiJ2,J 3 J4 „-i(J3||C (fe,B) ( m )ll J 4)n-i, (n > m) where the f n owing pr0 perty for the 6-j symbol with an 



element equal 1 has been used 



where (T k ) JlJ3> j 3 j 4 is defined in @. The term T^ 1_m_1 
in ( ftl| ) results from the iteration of (100) until one gets 
n = m. In the case when n = m in (|99|) we should apply 
( |9l| ) obtaining 





A 2 






/ A 2 


Ai 






J 2 


A] 




I J 2 


Ji 


i} 



^JslJ^O^in) ||(A 2 J 4 ),J 2 ), 



(101) 



^ (-1) A ^+ J »+ J »+* v/(2J 1 + l)(2J 2 + l) 

A ,i A J ;i 3 } n (Ai||^(n)||A 2 )„ 
Introducing ( |l0l|) into ( |98|) we get, 



APPENDIX C: THE STRING ORDER 
PARAMETER OF SPIN 1 CHAIN AND LADDER 

Let us first consider the spin 1 chain. The MP ansatz 
is given simply by, 



|JiM 1 ) JV =^ A JiJ* l(U 2 ), J1M1) 
J2 



N 



(107) 



„(Jiiic (fc ' B) wP2)n = E a 2£z 3 



^3 



where 



where the state |(1 J 2 ), J\M\)n reads as in ( |36| ) with 
(102) A = 1. We shall choose half-integer values of the an- 
gular momenta J± and J 2 which amounts to have a spin 
1/2 at one end of the chain H|t]|. 
We shall next show that the operators T = 1 and e l7rS * 



O 



(fc,B) _ ( A Xl 1 A x ' 2 (—i) x i+ J 2+J3+k Qnqhave both an eigenvalue equal to 1. Let us first of all write 

•^1^2,^3 ^AlA 2 V JlJ3 j ^2^3^ ' ^ ' . , , . 



x v /(2Ji + l)(2J 2 + l) j ^ * ^ 3 j (AiHO^HAa) 
This ends the proof of eq.(|4l|). 

Recursion relation of the energy expectation 
values 

We shall not give here the explicit proof of eq. (^0|) since 
it is quite analogous to the one performed in the previous 
paragraph. We shall simply state the result. 

The matrix hj u j a appearing in j50| ) is given by the 
sum 



out explicitely their components, 



(108) 



hj u j a = hy> t j 2 +h 



(2) 

Jl,J2 



(104) 



(T) 

Ji Mi J[ M[ , J 2 M 2 J 2 M 2 
— 8m 1 -M 2 ,M' 1 -M 2 ^JiJ 2 Aj' i j 2 
x (lMi - M 2 , J 2 M 2 \J 1 M 1 )(IM[ - M 2 , J^M^\J[M{) 

V / Ji Mi J[ M[, J 2 M 2 J' 2 

= 5 Mi -m 2 ,m[-m> {-l) M -- M * A, h . h Aj,j 2 (109) 
x (lMi - M 2 , J 2 M 2 | JiMi)(lAf( - M 2 , J^il^| J[M[) 

The normalization conditions on Aj ± j 2 read 
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J2 A 2 JxJ2 =l, VJa 
J 2 



(110) 



Let us come back to eq.(|ll5|), which can be written as 



Using these eqs. and the properties of the CG coeffi- 
cients, one can verify that v and v Bt defined as 



E 



-pjt 



a 



2Ji + 1 JlJ2 



(-1) 



Mi -1/2 



<(U 2 )J 1 M 1 |^|(lJ 2 )J 1 Af 1 ) 
(120) 



v JiM 1 J' 1 M[ — $JiJ[ &MiM[ 



(111) 



where Sf denotes the spin operator acting on the spin 1 . 
Using the Wigner-Eckart theorem we get, 



are right eigenvectors with eigenvalue 1 of the matrices 
T and e l7TSz respectively. 

Similarly the left eigenvectors associated to this eigen- 
value are given by, 

Pj 1 m 1 j[m[ = 5miM[ PjJ{2Ji + 1) (112) 
P S J lMl j>M> = h hM[ {-l) Ah - 1/2 pjJ{2Jx + 1) 

where pj is the left eigenvector with eigenvalue 1 of the 
matrix T JlJ2 = A^j^- 

According to eq. (E3) the string order parameter 17(00) 
is given by the product of two matrix elements which we 
compute below. 

Let us first consider, 

(113) 

= J2PJiMiJ' 1 M' 1 S z JlMl J[M[,J 2 M 2 J 2 M 2 v %m 2 j 2 m 2 

The hated version of S z is given by, 
(&) 

\ / Jj Mi J[ M[ , J 2 M 2 ,T 2 M' 2 

= S Mi -m 2 ,m[-m 2 {Mi - M 2 ) A JlJ2 Aj/j/ 
x (IMl - M 2 , J 2 M 2 | JiMi)(lM{ - M2, J 2 M^| J(M{ 

which together with ( 112| , 113| ) lead to, 



(-1) 



J 1 -Mi 



((U 2 )J 1 M 1 |5f|(lJ 2 )J 1 A/ 1 ) 

Ji 1 Ji 
-Mi Mi 

Mi 



x /Ji(2Ji + l)(Ji + l) 



((l^JiHSilKl^Ji) (121) 
((lJaJJxIISilKUaJJi) 



The reduced matrix element appearing in (121) can be 
computed using (|9l|), 



((l^JiHSiiKi^Ji) 

1 Ji j 2 
Ji 1 1 



= V6(-l) Jl+Ja (2Ji + l) 



(122) 



_ V2Ji + l(2+./i(.7i+l)-./ 2 (J 2 + l)) 
2 vOiCJi + l) 



Substituting (121 122) into (120) and performing the 
sum over Mi with the aid of the formula, 



E M (-!) 

M=-J 

(114) 

we get finally, 



Af-s 



(J + ~) (-1) J ~ 1/2 , (J : half integer) 



(123) 



(p\S*\if 



zE^ 



hj 2 



-i) 



J1-1/2 



(p\S>\^) = ErfSiA\j a (H5) 

X (_l)M 2 -l/2 (Mi _ Ma j (<1Mi _ J 2 M 2 | J lMl ))2 



2 + Ji(J! + l)- J 2 (J 2 + 1) 

JiOA + i) 



(124) 



Similarly we get 



From cqs(119 124) we immediately get the value of 



(p st \s z \v) = E 



P'l 



■A, 



2./1 + I ^Ji-h 

X (_l)Mi-l/2 (Mi _ Ma j (<1Mi _ J 2 M 2 | J lMl )) 2 

Observing that 

(_ 1 jJWi-l/a (Mi _ Ma j = _(_ 1 )M 2 -l/2 _ Ma ) 

(117) 

where Mi — M 2 = 0, ±1, we obtain 

( P st \S*\v) = -(p\r*\v st ) (118) 
which in turn implies 



g(oo) in the AKLT case, 
(116) 

AKLT : An = 1 -> .g(oo) = -(2/3) 2 



(125) 



In ref § the spin 1 Heisenberg chain was studied with 
a MP ansatz built up with two states with J = 1/2 and 
3/2. The values of the MP parameters obtained in Q] 
are reproduced below 



A11 = 0.988995, An = 0.14795 

2 2 2 2 

As 1 = -0.952887, A3 3 = -0.303325 

2 2 2 2 

pi = 0.97646, pi = 0.023539 



(126) 



ff (oo) = -((p|5'|« rt )) 1 



Introducing (126) into cqs.(119) and (124) we get 
(7(00) = —0.387, which can be compared with the ex- 
Ill 9 ) act value given by _o.374325 Q. In |] the spin 1 chain 
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was studied with a MP ansatz with two spin 1/2 and 
two spin 3/2 states, which yields 17(00) = —0.3759. This 
shows again that MP ansatzs with multiplicity improve 
considerably the accuracy of the numerical results • 

Let us go now to the spin 1 ladder. In section IV 
we gave an intuitive argument which suggested that the 
LRTO of the single spin 1 chains is destroyed by the 
interchain coupling. Next we show that this is indeed 
what happens. 

Let us first write eqs p8UlO^ ) in the case of ladders. 



(T) 



J 1 M 1 J' 1 M[,J 2 M 2 J 2 M 2 
X 



- S Ml -M 3 ,Mi-M^ Ea A JiJ 2 A J[J 2 ( 127 ) 

x (AMi - M 2 , J 2 M 2 \ JxMx){\M[ - M 2 , J 2 M 2 \J[M[) 



J 1 M 1 J[M[,J 2 M 2 J' 2 M 2 

Lav A \j 2 A*, j, (128) 
x (LMi - M 2 , J 2 M 2 \ JiMi){lM[ - M' 2 , J 2 M 2 \JiM{) 
x (AMi - M 2 |e i7rS J|A'MiM 2 ) 

where Sf denotes the spin operator acting on the first 
leg of the ladder. The vector vj 1 m 1 j , 1 m' 1 given in ( 111 ) is 
an eigenvector with eigenvalue 1 of the matrix T defined 



by (127). This property is a consequence of the normal- 
ization condition (39). For the spin 1 chain the latter 
condition also guarantees the existence of an eigenvalue 



1 of the operator (109). However this is not generally the 
case for the operator (128). 

can be deduced ex- 



The last matrix element in ( |125 
pressing the state |A//) of the rung in terms of the spin 1 
states of every site, 

|A/i) = \sm 2 ) 2 {lmilm 2 \Xn) (129) 

m\ra 2 

We thus get 



J 1 M 1 J[M[,J 2 M 2 J 2 M 2 
= ^Mi-M 2 ,M[-M' 2 J2xX' mi m 2 A J X J 2 A J[ J' 2 (~ l)™ 1 

x (IMl - M 2 , J 2 M 2 \.J 1 M l )(lM' l - M 2 , J 2 M 2 \J[M[) 
x (lmilm 2 |AMi - M 2 ) (lmilm 2 |A' 'M x - M 2 )) 



W e can actually set up M\ = M[ and M 2 — M 2 in 
( |130| ) since in the computation of the string order pa- 
rameter, the third component of the angular momenta is 
preserved. We have computed the highest eigenvalue x st 
of the matrix ( |130| ), which turns out to be smaller than 
one. This shows that g(£) decays exponentially with a 
correlation length £ st whose value is obtained by the eq. 



FIGURE CAPTIONS 

Fig. 1 Graphical representation of the MP ansatz ( |35| ) 
in the case of the spin 1/2 ladder and basis | JM)m with 
J = and 1. Every dot represents a spin 1/2. A link 
between two dots denotes the formation of a singlet be- 
tween the spins. Doted lines denote symetrization of the 
spins encircled by them. 

Fig.2 The MP parameters for the ladder AAi/ 2 - In 
figures 2 to 10 we adopt the notation x = |Jii/Jl|. The 
curve A j j 2 is labelled as [ J\ , J 2 , A] . 

Fig. 3 Same notations as in fig. 2 but for the FAi/ 2 
ladder. 

Fig. 4 J U (AF) = Jf{AF) and J U {FA) = J¥(FA). 

Fig. 5 GS energy per site of the AA and AF ladders 
and their U dual models. 

Fig. 6 Spin-correlation lengths of the A A ladder and 
its U dual. 

Fig. 7 GS energy per site of the AF ladder and its T 
dual given by the FA ladder. 

Fig. 8 GS energy per site of the AA ladder and its S 
dual. 

Fig. 9 Spin-correlation length of the AA ladder and its 
S dual. 

Fig. 10 Plots of the spin-correlation length £ and the 
string correlation length £ st of the ladder AA\. 

Fig. 11 Pictorical representation of a possible AKLT 
state of the spin 1 ladder. 

Fig. 12 Graphical representation of a generic AKLT 
state of a ladder denoted as (p, q, r). There are a total of 
p + q + r dots inside every circle representing a total spin 
S = (p + q + r)/2. 



(130) 



f* = -l/ln\x st \ 



(131) 
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